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^ ! Abstract 

f~^ i The active-sterile neutrino conversion probability is calculated for neutrino 

^^ , propagating in a medium in the presence of random magnetic field fluctuations. 



Necessary condition for the probability to be positive definite is obtained. Using 
_ . this necessary condition we put constraint on the neutrino magnetic moment 

D , from active-sterile electron neutrino conversion in the early universe hot plasma 



p. 



X 



and in supernova. 
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1 Introduction 

Neutrino propagation in the presence of a strong magnetic field is of great interest 
from the point of view of astrophysics as well as early universe hot plasma[|l|. The 
presence of primordial magnetic fields {B c:^ lO^^G) over galactic scales, when ex- 
trapolated back can give very large fields and these large fields might have affected 
the particle interactions in the early universe0]. It is believed that the primordial 
plasma was chaotic in nature and the magnetic flux lines moving along with this 
plasma might get mixed up creating randomness in the fields itself. It is also be- 
lieved that the magnetic field inside a newly born neutron star is huge and random 
in nature p, ^, ^, Q. This strong random magnetic fields can influence the neutrino 
propagation as well as their conversion in the early universe hot plasma and in a 
supernova mediump], ^ |^, ^. Neutrinos having magnetic moment or transition mag- 
netic moment can flip their chirality and left handed neutrinos becomes right handed 
sterile neutrinos in the presence of an external magnetic field. These right handed 
neutrinos being sterile to weak interaction stream out from the sun or supernova core. 
The mechanism of helicity flipping in the presence of a magnetic field is one of the 



explanation for the supernova cooling mechanism or the solar neutrino deficit [jTO|. 
Recently Pastor et. al.^] have derived bounds on the transition magnetic moment 
of the Majorana neutrinos in the presence of random magnetic fields from the super- 
nova energy loss as well as from nucleosynthesis. In this paper we discuss about the 
active-sterile electron neutrino conversion {u^l -^ ^eit) in the presence of a random 
magnetic field in early universe hot plasma and in supernova medium. From the 
positive definiteness of the average conversion probability we put constraint on the 
neutrino magnetic moment. 

The paper is organized as follows: In Section 2, we have derived the average 
probability equation for the active sterile/active {va ~^ ^x) neutrino conversion in 
the presence of a random fluctuation over the constant background magnetic field. 
For this derivation we have used the simple delta correlation for the uncorrelated 
random magnetic field domains. We show that, the fluctuation in the transverse and 
longitudinal components of the magnetic field are mixed up. We obtain the solution 
for the probability equation and found the necessary condition for the conversion 
probability to be positive definite. We have considered the neutrino propagation in 
the early universe plasma and in the supernova medium. Assuming the magnetic 
fields in the early universe hot plasma and inside the core of a newly born neutron 
star to be strong and purely random in nature we put constraints on the magnetic 
moment for active sterile neutrino conversion, which are discussed in Section 3 and 
4 respectively. In discussion Section 5, we briefly summarise our results and discuss 




about the further improvement of the bound on the neutrino magnetic moment. 

2 Neutrino propagation 

The wave equation for the propagation of two neutrinos, one active and one sterile 
(or active) in a plasma in the presence of a magnetic field is given by 

where x = s,b (sterile, active) and 6 = e,/i and r. In the evolution eqn.(|l]) the 
diagonal components are Haa = Kec — A + Vaxiai and H^x = and the off-diagonal 
entries Hax = Hxa = fJ'B±{t). The A is given as A = cos(2^) {1712 — m\)/2E and 
vanishes for degenerate Majorana neutrino and E is the neutrino energy. For Dirac 
neutrino A = and /x corresponds to its diagonal magnetic moment. On the other 
hand for Majorana neutrino /i is the transition magnetic moment. V^ec is the differ- 
ence of neutrino vector interaction potentials, Kec = Kec ~ Kec- Here we consider 
neutrinos without mixing, which implies cos(26') = 1. The axial vector potential 
Vaxiai = fJ'eff^-^/k is generated by the mean axial vector current of charged leptons 
in an external magnetic field 0, |Tl], |12[ . 

Let us define the functions R = Re{< v^Vx >) and / = Im{< v^Vx >). Then 
using these in eqn. (|l|) we obtain the following set of equations: 

dP{t) 



dt 

dl{t) 



-2Hax{t)I{t), (2) 



Haait)R{t) + Hax{t){2P{t)-l), (3) 

dR{t) 



dt 
and 

-Haa{t)I{t). (4) 



dt 

The function P{t) is the neutrino conversion probability Pu^^i,^(t). The magnetic 
field can be written as B{t) = Bq + B{t), where Bq is the constant background field 
and B{t) is the random fluctuation over it. Putting this in Haa{t) and Hax{t) we get 

Haa{t) = HaaiO) + H,a{t) = (Kec " A + /ie//fi||o) + f^effB\\{t) (5) 

and 

Hax{t) = HaxiO) + Hax{t) = fiB^o + f^B^{t). (6) 

For the neutrino conversion length greater than the domain size [Iconv >> -^o) 
one can average the the equations(H), (^ and (^ over the random magnetic field 



distribution. Let us define the average of tlie functions < P{t) > = V(t), < 
R{t) > = TZit) and < I{t) > = T{t). Using these in eqns.(|D and (H) we obtain 

^ = <if,.(t)(2P(t)-l)> + <iJ,„(t)i?(t)> 
at 

= HaM < (2P(t) -!)> + < Ha.{2P{t) - 1) > 

+Haa{Q) < R{t) > + < Haa{t)R{t) >, (7) 

and 

^^ = -Haa{Q) < lit) >-< HUt)m > (8) 

respectively. The magnetic field in different domains is randomly oriented with re- 
spect to the neutrino propagation direction. So the neutrino conversion probability 
depends on the root mean square [rms) value of the random magnetic field. With 
the use of the delta correlation for uncorrelated magnetic field domain of size Lq, the 



average of the random magnetic field is|13, 14, 15 



< B\\{t) >=< B±{t) >=< B\\{t)B±{t) >= 0, (9) 

< 5,||(t)5,||(ti) >=< Bf > 5ijLo5{t - ti), (10) 

and 

< 5ix(t)5,±(ti) >=< Bl > S,,LoS{t - ti). (11) 



The rms value of the averaged magnetic field is given as Brms = V< B"^ >. Neglecting 
the higher order correlation of Haaif) and Haxif) with /(t), R{t) and P{t) and using 
the above equations in eqn.(|^) to (|) we obtain 

^ = Haa{Q)n{t) + HaMi'^nt) " 1) " 2(1^ + ri,)J(t), (12) 

dR{t) 

dt 

and 



-if„,(o)J(t)-2^||7^(t), (13) 



dt 
For convenience we define 



and 



^^(^) ^ -2/7,,(0)X(t) - r^ (2P(t) - 1) . (14) 



where 



X(t) = e-2(r[i+rx)tx^(t)^ (15) 

7^(t) = e-2^ll*7^l(t), (16) 

r± = ^fi' <B^> Lo, (17) 



and 



1 



II = e^e// 



2 <B^ > La. 



r± and rii are the transverse and longitudinal magnetic field damping parameters. 



Putting these in eqn.(14) and differentiating two times with respect to t we obtain 



d^Vjt) 



+ 4(rx + r 
+ 



d^Vjt) 

dt^ 



+ 4 sr^rii + ri + r^ + h^,{o) + 



rir,, + T^Tl + H^M^w + ^" , ^ ) v{t) 



4 rir|| + r^r^ + if2^(o)r|| + 



HLm_ 



0, 



H!M \ dv{t) 

4 / dt 



(19) 



dV{t) I 



r_L and ^-^/i=o 



with the boundary conditions V{t)/t=o = 0, 
2H^^{0) — 2T\. Switching off the damping terms in eqn.(|^), will give MSW type 
solution with magnetic field[0, |10|, |l|]. For Haa{0) = Hax{0) = T\\ = 0, the equ.® 
reduces to 

0, (20) 



^^^W+4r,^ + 4riP(t)-2rl 



dt^ ' "^ dt 
and it has the same solution as shown by Pastor et. al.|P]. In this the effect of strong 

random magnetic field on the neutrino transition magnetic moment is studied both 



in the early universe hot plasma and in supernova. The eqn.(|T^) can be written in a 
simplified form as 



'^^^)^A.m,B/-m^2c.rit)-c.^o, 



dt^ " dt^ 

where the quantities Aq, Bq and Cq are 



dt 



A, 



4 (Fx + F| 



Bo = 4 3FxF|| +TI + TI + HIM + 



HIM' 



and 



Co = 4 FiF|| + F^F^ + HlM^W + 



HLm. 



(21) 

(22) 
(23) 

(24) 



respectively. Eqns.(^) and ( p4D shows that the damping terms are mixed up among 
themselves and also with the terms -ffaa(O) and HaxiS^)- The solution of eqn.(^) is 

1 



where y{t) is given as 

y{t) = 






Vit) = y{t) + 



Ai (e^^i* - cos(Z4t) - ^^ sin(Z4t) 

I 2 Z 4 



(25) 



cos(Z4t) , (F^-f -^ 



3Zi 

Z4tj - 
-) 

^— sin(Z4t) 



(26) 



The coefficient Ai is given as 



A, 



{Al - UAqT^ + 36ri - SGHlJO) + SApZ^ - IST^Zi + \Zl + QZ; 

{^Zl + 18Z|) 



and Zi and Z4 are 
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'2 



V 2 



+ 



(f) 



f)^(?)T 



(27) 



(2S 



and 



/3 






(!) 


2 


- V 


2 


1 


W(ir+«T 



1 

2 



respectively. The quantities p and g are defined as 



P 



42 



and 



g 



(2Co-^+ ^ 



Ai 



py (q 
3) 



(29) 



(30) 



(31) 
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The solution eqn.(^) for the average neutrino conversion probability V{t) is very 
complicated. So here we will consider condition for the existence of the solution 
rather than going into the details of it. Since the probability is positive definite 
(0 < V(t) < 1), the following condition, 



h 



py fq 

3) \2 



> 0, 



(32) 



has to satisfy, otherwise Zi and Z4 will be complex so also V{t). Then putting values 
of p and q in eqn. (p2|) we obtain the condition 



^Hlio) + Hl{o) > l(ri + rj-r^r||). 



(33) 



Thus for the average neutrino conversion probability to be positive definite, the above 
condition has to satisfy, irrespective of the form of the neutrino potential and the 
magnetic field. Let us consider the neutrino propagation in the early universe hot 
plasma and the core of the newly formed neutron star, where the above condition 
can be satisfied. 



3 Early Universe hot plasma 

In the early universe when the temperature T » 1 MeV, all the particles are in 
thermal equilibrium [T^. In the early universe the magnetic flux lines are moving 
along with the hot plasma. Because of the chaotic motion of the plasma, the flux 
lines will be mixed up and twisted thus creating the randomness in the magnetic 
field. So the constant background field Bq must be very small compared to the 
random fluctuation term. Thus we can safely assume i?o — and then Hax{0) — 
and fieffB\\Q ~ 0. For small particle-anti particle asymmetry in the early universe hot 
plasma the axial vector potential contribution can be very small as in a relativistic 
plasma the charged lepton/anti-lepton masses and their chemical potentials are small 
compared to the kinetic energy term, thus the factor fieff will also be small^, 0]. So 
we can neglect the longitudinal damping term for the magnetic field. Then eqn.(^3D 
will be 

HUO) > ^r^. (34) 



The vector potential for a neutrino in the early universe hot plasma isfTg 

T2 



Kec = V2Gpn^{T) 



L-A 



Mir, 



(35) 



Here n^{T) ^ 0.244T^ is the photon number density and A ^ 55. For temperature 
T >> rrie the second non-local term is greater than the first term for very small 
particle-anti particle asymmetry and the second term is 
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Kec ^ -3.45 X 10-20 (^^^^-^) MeV, (36) 

for electron neutrino. For z/g^ -^ UeR process we assume that the right handed 
neutrino produced is sterile and decouple from the system. 

Irrespective of its origin, the primordial magnetic field had a very large value 
in the early universe. We assume that the primordial plasma consists of magnetic 
domain structure with a size Lq and the magnetic field is uniform and constant 
within each domain and the field in different domains are randomly aligned H, ITB 



For homogeneous magnetic fields, the flux conservation indicates that B oc T^. But 
detailed structure of random magnetic field profile in the early universe hot plasma 
will depend on the complicated nature of the magneto- hydrodynamic equations ^^ 



For the root mean square field Brms = V< B"^ >, averaged over a volume L^ >> L\ 



3 



we assume a power law behavior 1 20, 21 



Brms = Bn[—\ ( — j , (37) 



where Tq is the temperature at some reference epoch and Bn is the corresponding 
field strength within a domain of size Lq. The maximal scale we have chosen is the 
horizon length L = Ih = Mpi/T'^, where Mpi = 1.22 x 10^^ GeV is the Planck mass. 
We take here BnT^L^ to be constant. If the primordial magnetic fields are to be 
the seed fields for the galactic dynamo, then it should survive till the recombination 



epoch and this leads to a minimal domain size||22||, 

Lo > L™" ^ 10=^ cm X (^^) . (38) 

Here we take i?„ = 10^^ Gauss and Tq = Tew = 10^ MeV the electroweak phase 



transition temperature 0]. Putting all these in eqn.(0) we obtain 



1.15 X 10-15+^2n . J, (l^„)^ (i+„) 

The index parameter n = corresponds to the uniform magnetic field, which is not 
physically very likely in the early universe hot plasma. For the random magnetic field 
along the neutrino trajectory the index parameter would he n = 1/2 and n = 3/2 
corresponds to 3-dimensional elementary cells[^, |], [1^. As the magnetic field is 
random in nature we take n = 1/2. For QCD phase transition temperature T = 
Tqcd ^ 200 MeV we obtain /i < 8 x IO-^Vb- 

4 Supernova Core 

Now let us consider the neutrino propagation in the supernova medium. As shown 
in eqn.(|D and eqn.(|) we have -ffaa(O) = (Kec - A + ^effB\\o) and -ffax(O) = /i5_Lo- 
Thomson and Dunkan]^ have shown that, magnetic fields as strong as 10^^ to 10^^ 
Gauss might be generated inside the core of the supernova due to a small scale dynamo 
mechanism. If these fields are generated after core collapse, then it could be viewed 
as random superposition of many small dipole of size Lq ~ 1 Km. So we neglect 
here the Bq term in the diagonal and non-diagonal parts. Thus Haa{0) = (Kec — A) 
and Hax{0) = 0. For u^l -^ ^eR the vector potential experience by neutrino in the 
supernova medium is 

Kec = 4 X 10-'puf{Y,) MeV, (40) 

where f{Ye) = {3Ye — 1) is the electron neutrino abundance factor and pi4 is the 
density in units of 10^^ g/cm?. The right-handed electron neutrino being sterile 
stream away from the supernova core. The A = 5 x 10~"^^/ii^ioo (Am^/e\^^) MeV, 
where the neutrino energy Eiqq is in units of 100 MeV. We can see that A << Kec 

pvpn fnr larp"p Am TnsiHp thfi npiitrnn star rnrp thprp arp Ipss nnmbpr nf plfiptrnns 



so the damping term Fy can be very small compared to the vector potential or the 
transverse damping term F^. Assuming Fy to be very small inside the core we have 
Haa{0) ^ Kec- Then eqn.(^3|) is reduced to the same condition as in the early universe 
hot plasma -ffaa(O) > 2T±/^/3. Inside the core the damping term is 

F, = 2.2 X 10-Vo^ (%^) (^) MeV, (41) 

where /io is in units of fj,B the Bohr magneton. Putting the values of Haa{0) and F^ 
in eqn.(0) we obtain for the magnetic moment 

( Brms \ I i-IQ \ 

\ G ) \cm) 

For Brms ~ 10^^ Gauss, Lq ~ 1 Km, pi4 ~ 8 and the electron abundance factor 
Ye ~ 0.3 inside the neutron starQ, the magnetic moment comes out to be yU < 
3.5 X IQ^^^iiB- For smaller value of the magnetic field, the magnetic moment will be 
larger. 



5 Discussion 

Assuming the magnetic field has a random fluctuation over the mean value, we have 
derived the average probability equation for neutrino conversion/spin precession in 
the magnetized plasma. As a consequence of the random fluctuation in the magnetic 
field the transverse and longitudinal magnetic field damping are getting mixed up. We 
have assume a delta correlation for the random magnetic field domains and obtain the 
solution for the probability equation. The definiteness of the probability is invalid 
if the condition in eqn.(p3D is not satisfied. We assume the magnetic field to be 
purely random in nature, in the early universe hot plasma as well as in supernova 
medium (inside the newly born neutron star core) and consider the process v^l ~^ '^eR- 
At the QCD phase transition temperature T ~ 200 MeV we obtain the neutrino 
magnetic moment /i < 8 x 10~^^ fiB- Inside the supernova core we consider the same 
process UeL — > ^eR and for Brms — 10^^ Gauss, obtain /x < 3.5 x lO^^^/i^. In this 
estimate we have neglected the contribution from the axial vector potential as well 
as the contribution due to the constant magnetic field in the plasma. We have also 
neglected the contribution from the longitudinal damping term. So inclusion of all 
these terms might improve the constraints on the magnetic moment. Apart from 
that the magnetic field profiles for the supernova and early universe are very much 
speculative, so the upper limit might change for other magnetic field profiles. 
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